Stability of topologically invariant order parameters at finite temperature 
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Topological entanglement entropy is a topological invariant which can detect topological order 
of quantum many-body ground state. We assume an existence of such order parameter at finite 
temperature which is invariant under smooth deformation of the subsystems, and study its stability 
under hamiltonian perturbation. We apply this assumption to a Gibbs state of hamiltonian which 
satisfies so called 'strong commuting' condition, which we shall define in the paper. Interesting 
models in this category include local hamiltonian models based on quantum error correcting code. 
We prove a stability of such topologically invariant order parameter against arbitrary perturbation 
which can be expressed as a sum of geometrically local bounded-norm terms. The first order 
correction against such perturbation vanishes in the thermodynamic limit. 



In a quantum many-body system at zero temperature, 
system can exhibit topological order depending on the 
parameters of the hamiltonian. Prototypical examples 
include Kitaev's quantum double model and Levin- Wen 
string-net model^^ Unfortunately the interaction terms 
of these models are usually more than two body. This 
might be problematic, since it is hard to find materials 
which can be described exactly by such models. One may 
hope, however, that these hamiltonian arises as an low- 
energy effective hamiltonian of a simpler model. If the 
perturbative terms in the effective hamiltonian is small 
enough, one can argue that the topological order may 
still be protected^ - — For gapped frustration-free hamil- 
tonian, certain topological order condition on the ground 
state ensures a gap protection against geometrically local 
bounded norm perturbation. 

It is plausible to conclude gap protection implies the 
stability of topological order, since absence of quantum 
phase transition is likely to be an evidence that they be- 
long to the same phase. This logic cannot be applied to 
finite temperature, since the density operator has non- 
trivial support on excited states as well. This problem 
is not so interesting in 2D, since there is no finite tem- 
perature topological order in 2D£~— However, nontrivial 
topological order might survive at finite temperature for 
higher dimensional system. For instance, 3D generaliza- 
tion of Levin- Wen type topological entanglement entropy 
is nonzero below the critical temperature of 3D toric 
code£ Other models with similar structure will likely 
show a similar behavior This type of topological or- 
der is of classical origin. Even though the ground state 
of the system encodes nontrivial quantum information, 
only classical information survives at finite temperature. 
4D toric code on the other hand can successfully preserve 
quantum information even at finite temperature^ 

One of the main objectives of this paper is to estab- 
lish an analogous stability result at finite temperature, 
potentially to the systems discussed above. We have ap- 
proached this problem by observing the common proper- 
ties of these systems. When the system is topologically 
ordered, there is a nontrivial topological entanglement 
entropy. Important property of topological entanglement 



entropy is that it is invariant under small deformation of 
the subsystems as long as they do not change the 'shape' 
of the configuration. When the system is in a topologi- 
cally disordered state, topological entanglement entropy 
is likely to be 0, as evidenced in the calculation of 3D toric 
code<2 The topological entanglement entropy of this sys- 
tem is still 'invariant' in a sense that it stays with a 
small deformation of the subsystem. In an actual model, 
such invariance is more likely to be an approximation of a 
small number that vanishes in the thermodynamic limit. 

The paper can be roughly divided into two parts. First 
part is about studying the implications of such topolog- 
ical invariance. We shall present a general argument 
that such systems must not have any long range or- 
der. Furthermore, we shall show that such invariance 
implies an asymptotic conditional independence for cer- 
tain configurations. A tripartite state pabc is condi- 
tionally independent if a conditional mutual information 
I(A : B\C) is 0. Asymptotic conditional independence 
means I (A : B\C) = e for some e that vanishes in 
the thermodynamic limit. It is important to note that 
asymptotic conditional independence does not hold for 
arbitrary configurations. We shall study the configura- 
tions arising from the invariance of topological entangle- 
ment entropy under smooth deformation. 

The rest of the paper is about showing the implication 
of asymptotic conditional independence. We shall define 
a class of models which satisfy the condition called strong 
commuting condition. For such models, one can bound 
a first-order perturbation of Levin- Wen type topological 
entanglement entropy under a sum of geometrically local 
finite-norm perturbations. The implication of this result 
is twofold. If there exists a topological order in a sense 
of having nonzero topological entanglement entropy, it is 
stable against hamiltonian perturbation. Hence the topo- 
logical state of matter is robust. On the other hand, it is 
impossible to create topological order from a disordered 
system by adding a small perturbation. 

The proof depends on two statements about quantum 
many-body systems. The first statement is a certain vari- 
ant of Lieb-Robinson bound. This technique was used 
by Hastings in proving certain locality properties of fi- 
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nite temperature quantum systems i 17 ' 18 Second state- 
ment concerns a spectrum of an operator whose trace 
reduces to conditional mutual information. The second 
statement is the only part where the strong commuting 
condition is needed. Logic of the rest of the paper re- 
mains intact as long as the hamiltonian can be expressed 
as a sum of geometrically local bounded-norm terms. We 
shall discuss how the second statement may potentially 
be generalized to a wider class of models. 

Once these technical results are established, the con- 
ceptual idea behind the proof is quite simple. Effect 
of the perturbation by local terms can be written as a 
sum of various correlation functions. These terms can be 
bounded by conditional mutual information, which was 
assumed to be small. We set the stage by motivating 
and precisely formulating the asymptotic conditional in- 
dependence condition for topologically ordered system in 
Section |TJ We sketch proof and present the mathematical 
ideas in Section [TTJ The main stability result is stated in 
IIIII Section ITVl will discuss how this result can be applied 
in the contexts of known models. 



I. ENTANGLEMENT ENTROPY OF 
TOPOLOGICALLY ORDERED SYSTEM 

Ground state of topologically ordered system satisfies 
an area law. 

S A = c\dA\ -7 + o(l), (1) 

where 7 is a quantum dimension of the system. It is 
usually implicitly assumed that A is a simply connected 
subsystem of the entire lattice A. If A changes into a 
topologically distinct object, 7 should change as well. 
For instance, if A is a union of two simply connected 
subsystems which are widely separated, 7 should be re- 
placed by 27. This tells us that 7 must be not only a 
function of the state, but also a function of the topol- 
ogy of A, hence we write it instead as ja from now on. 
Unfortunately topology in the spirit of classifying shapes 
of smooth objects are not a well defined concept for dis- 
crete systems. For instance, consider a torus represented 
as n x n grid, both ends identified. Open string can be 
mapped into a closed string and vice versa via 'smooth 
deformation,' which corresponds to adding or subtract- 
ing a single cell. One remedy for such problem is to 
embed the entire system onto a manifold and designate 
a finite neighborhood near each particles, so that the fi- 
nite neighborhoods partition the entire manifold. Under 
this embedding, one can easily see that open string and 
closed string have clear topological meaning: 'smooth de- 
formation' once lifted to the entire manifold cannot map 
open string into closed string and vice versa. This is the 
setting in which we shall explain our result. When we 
say subsystem A, it shall mean a subset of the manifold 
which is a union of finite neighborhoods near each parti- 
cle contained in A. Similarly, when we say an area and 



volume of A, it shall mean the area and volume of the 
embedding on the manifold. 

If only depends on the topology of A, it must be 
equal to 74/ for two topologically equivalent subsystems 
A ~ A . The easiest nontrivial implication can be seen 
in FIGfTJ For each configurations, 74 is either 7 or 27, 
depending on if it is a simply connected region or union 
of two such widely separated subsystems. For all of these 
configurations, I(A : B\C) = o(l). 



A A A A 



B B B B 



(a) (b) (c) (d) 

FIG. 1: Three possible configurations of A, B and C that 
produces I (A : B\C) = o(l). When the subsystems are con- 
tingent to each other, they share a boundary. Otherwise, they 
are assumed to be far away from each other. The last figure 
denotes (7 = 0. 

In fact, knowing the exact value of 7^ is unimportant 
for showing such cancelations as long as we can find a 
pair of shapes which are topologically equivalent to one 
another. For instance, consider FIGHJ The constant sub- 
correction term of Sab and Sb are both 7, so they cancel 
out each other trivially. ABC and BC are both a closed 
surface with a puncture. Since their shapes are topolog- 
ically equivalent, the constant subcorrection term cancel 
out. One can easily check that the area term cancels out 
each other as well. 

Equation[T]is a strong evidence that degrees of freedom 
is localized near the boundary of A. In general, what is 
more likely to happen is that the majority of the degrees 
of freedom are located near the boundary and there are 
small tails that decay as we move away. Due to this de- 
cay, tails are negligible when A and C are far away from 



A 



C 
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FIG. 2: Another configuration which yields I(A : B\C) = 
o(l). 
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each other, but it may not be when they are close to- 
gether. From these observations, we arrive at the precise 
definition of asymptotic conditional independence. 

Assumption 1. If i)AB ~ B and BC ~ ABC or 

ii)AB ~ ABC and BC ~ B, I (A : B\C) < f{l AC ) 
for some decaying function f , where Iab is the distance 
between A and B . 

Volume term emerges at finite temperature, but the 
cancelation property remains the same: Vol(AB) + 
Vol{BC) - Vol(ABC) - V{C) = when A,B and C 
only coincide on the boundary. Assumption [1] implicitly 
assumes a finite correlation length. Hence we do not ex- 
pect it to hold for critical systems. 

If the configurations for extracting the topological en- 
tanglement entropy is known, there is an alternative way 
to see the asymptotic conditional independence. Con- 
sider a Levin- Wen type configuration or its 3D general- 
ization. Topological entanglement entropy can be writ- 
ten as I (A : B\C) = 2j. If this quantity is a topologi- 
cal invariant, it must remain the same under a deforma- 
tion that preserves the shape of the subsystems. Sim- 
ple algebraic manipulation shows I (AD : B\C) — I (A : 
B\C) = I(D : B\AC). We have a small number on the 
left hand side, since the assumption was the invariance 
of topological entropy under small deformations. This 
may not be true when D is close to B. Hence right hand 
side must be a function that primarily depends on the 
distance between D and B which converges to 0. No- 
tice that the configurations arising from I(D : B\AC) 
in FIGJ3] are topologically equivalent to configurations 
arising from I(A : B\C) in FIGf2J Since 3D generaliza- 
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FIG. 3: Levin- Wen construction for topological entanglement 
entropy and its deformation. 

tion of Levin- Wen configuration can be also written as 
I(A : B\C) for some configurations, same logic can be 
applied as well. Notice that we did not assume anything 
about the specifics of the UV-divergent terms of entan- 
glement entropy. If the topological entanglement entropy 
is invariant under smooth deformation of the subsystems, 
conditional mutual information must be a small number 
in the thermodynamic limit. 



These two different perspectives have their own caveat. 
If we assume certain cancelation property of UV- 
divergent terms, the asymptotic conditional indepen- 
dence condition is obtained. On the other hand, as- 
suming the existence of certain topologically invariant or- 
der parameter automatically implies the same result irrc- 
gardless of the specifics about the entanglement entropy. 
It is not entirely clear, however, that if there always ex- 
ists such order parameter. Furthermore, the smallness 
of I (A : B\C) for configuration in FIG 1(d) is essential 



for our proof. One cannot derive it solely based on the 
topological invariance of certain order parameter. This is 
still a plausible assumption nonetheless, since the system 
likely possesses a finite correlation length in noncritical 
systems. 

Either way, these evidences strongly support the valid- 
ity of asymptotic conditional independence for topologi- 
cally ordered system. Remainder of this paper is about 
studying the consequence of this statement, but we need 
to be specific. We are trying to bound a first-order per- 
turbation of topological entanglement entropy under a 
sum of geometrically local finite-norm terms, but what 
kind of topological entanglement entropy should we use? 
Levin- Wen type configuration and its generalization seem 
to work, but the proof is rather general. For this reason, 
we would like to start with an abstract definition, and ob- 
serve how it can be applied to the configurations studied 
in the literature. 

We shall name linear combination of entropies in gen- 
eral as entropic order parameter S({Ai\) = ^aiSAt- 
Spatial deformation of the order parameter can be writ- 
ten as S({Ai}) — <S({A-}) in general for some deformed 
subsystem A^s. We are interested in a following class of 
entropic order parameters. 

Definition 1. An entropic order parameter S({Ai}) is 
(h,l2)-deformable if for all local site B, there exists de- 
formation {Ai} — > {A^} such that 



S({A l })-S({A'}) = J2*jI(3i 



(2) 



for dist(Ji, j'2) > W> mindist(A / i , B) > l^. {aj} is some 
sequence of numbers. 

Important property of the order parameters satisfying 
this definition is that one can always deform away the or- 
der parameter from any local site with distance at least 
lx- As we shall sec, effect of perturbation can be written 
as a sum of thermal correlation function between two ob- 
scrvablcs. For large enough distance, thermal correlation 
between the local site and the deformed order parameter 
becomes negligible due to the existence of finite correla- 
tion length. What remains to be seen is the thermal cor- 
relation between the difference between these two order 
parameters and the local site. There is no guarantee that 
the distance between those two observables to be large 
in general, but we shall show that it is possible to bound 
it nonetheless from conditional mutual information. The 
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perturbation bound vanishes in the limit li,fa — > oo, so 
one must check if the order parameter is Z 1; Z 2 -deformable 
for large enough value of l\ and fa- 



A. 2D and 3D Levin- Wen construction 

First example of entropic order parameter which allows 
a large value of fa and fa is Levin- Wen configuration in 
2D. Assume A and B are 3R x R rectangle, and C a 
union of two RxR squares. We can show that topological 
entanglement entropy I(A : B\C) is (-j, -j)-deformable. 
Let s denote a local site. If s is distance -j or larger away 
from ABC, the conditions are automatically satisfied. 
Let the distance from s to A, B, C be cIa, ds,dc- At least 
one of <1a or ds must be larger than ^. Without loss of 
generality, assume dA < ds- First consider deforming 
A into A\D. Using I(A : B\C) - I(A\D : B\C) = 
I(D : B\AC \ D). For any point, it is possible to choose 
D such that the distance between s and A \ D is larger 
than -j. For the distance between C and the s, the effect 
of deformation can be written as I(D : B\AC \ D) — 
I(D : B\C \D). Again, it is possible to choose D such 
that the distance between s and C \ D is larger than -j . 
An example of such deformation is illustrated in FIGfJ] 



— ( ■ 


A 




C 


B 



(a) 



(b) 



FIG. 4: Configurations before and after the deformation. 
Black colored region denotes the position of the local oper- 
ator. Irregardless of the position of this operator, one can al- 
ways deform the original operator so that it can be seaparted 
from the local operator. 

Similar idea can be applied to 3D generalization of Levin- 
Wen construction. See FIGJSJ As in the 2D case, one 
can deform away the region around a local site. We did 
not denote the position of the local operator due to the 
clarity of the picture, but it is located near the center 
of cubic deformation shown in FIG 5(c)| 5(f) 5(g)[ 5(1: 



There is an alternative 3D generalization for Levin- Wen 
model. This can be easily reconstructed by stacking slices 
of FIG|4]together. Same logic can be applied here as well. 



B. Alternative configurations 

Recently Grover et al. introduced a number of ways to 
extract the constant subcorrection term of entanglement 
entropy in 3DJ£ These are generally a 3D generalization 




(e)A'C' 



(f)B'C' 



(h)A'B'C' 



FIG. 5: One possible generalization of Levin- Wen construc- 
tion in 3D. Topological entanglement entropy 7 can be ex- 
tracted by I(A : C\B) = 27. The first four configurations 
denote the original configuration. The next four configura- 
tions denote the configurations after the deformation. 



of Kitacv-Prcskill construction in some ways. This means 
the topological entanglement entropy is extracted by a 
formula S top o — Sa + Sb + Sc - Sab — Sbc — Sac + 
Sabc f° r some configuration A, B, and C. One of their 
proposal, which is a three equipartition of torus allows a 
(h, y-deformability for large values of l\ and fa- 




FIG. 6: Alternative 3D configuration that allows deformabil- 
ity for large values of h,h- 



C. Some configurations are not deformable 

In showing the deformability of the configuration, it 
was vital to assume that at most two subsystems involved 
in the computation of entropic order parameter must be 
contingent on each local sites. There are several con- 
figurations which are known in the literature to disobey 
this rule. Primary example is Kitaev-Preskill construc- 
tion in 2D: there exists a triple point where three of the 
subsystems involved in the computation of entropic order 
parameter meet together. One of Grover et al.'s construc- 
tion also share a similar property. These are depicted in 
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FIG 13 These are naturally excluded from the scope of 
the stability proof. 
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(b) 



FIG. 7: Configurations which are not deformable with large 
values of h, fo- 



il. TECHNICAL RESULTS 

There are two main technical results that deserve to be 
discussed separately. First is a variant of Lieb-Robinson 
bound. Consider an infinitesimal perturbation V on the 
hamiltonian H , H —> H + eV. Such perturbation can be 
characterized by a directional derivative of H. 



Definition 2. 

d V H f{H) 



lim 

e-s-0 



f(H + eV)-f(H) 



(3) 



If the perturbation V is local, Hastings showed that 
1 



die-** 



-MV'e 



-I3H 



V), 



(4) 



for some quasi- local operator I/ VV 8 In the eigenbasis of 
H , V can be expressed as the following. 



V 



Tr tanhxi,- 
Va - 



(5) 



with X; 



-EA 



t ,j — 2 • Ei s are the eigenvalues of H. Im- 
mediate application of this formula concerns an effect of 
local perturbation on the thermal expectation value of 
some observable. 



d%(a)=R B (C(V' i a)), 



(6) 



where C(Oi, O2) = (O1O2) — (Oi) {0%) is a connected cor- 
relation function. One of the general recurring themes 
throughout this paper is that V and a are observablcs 
which are far away from each other. In those cases, it 
is possible to bound the connected correlation function 
from mutual information. Locality of V can be trivially 
shown for hamiltonian which consist of geometrically lo- 
cal commuting terms, but it is also important to note 
that its operator norm is bounded by ||V|| < \\V\\. This 
is explained in Appendix lAl In this context, locality of 
V can be stated in a following way. 



Lemma 1. Let r be the maximal radius of interaction 
of the local terms in H . If the local terms commute with 
each other, 



[V',O} = 



(7) 



for dist(supp(0) , supp(V)) > r. supp(0) is the minimal 
nontrivial support of O. 

The second idea concerns a special property of condi- 
tional mutual information. Note that conditional mutual 
information I (A : B\C) can be written in a following 
form. 



I(A:B\C)=tr[pH A:Blc ], 



(8) 



with H A , B \ C = In p abc + lnp B - hip^s - lnpsc- bip.4 
actually means ln(pA ® Ia^ ) — In Pa ® Ia c ■ Consider an 
observable a whose support lies on BC. For a classical 
system, 



\Re(tr[po-H A:Blc ]\<I(A:B\C)\\o-\\. 



(9) 



One way to see this is to take a partial trace over subsys- 
tem A. a is unaffected since it does not have nontrivial 
support on A. Using the fact that pabc = Pa\bcPbc 
and the positivity of relative entropy D(pa\bc\\pa\b), 
one can show that trA[pABcHA-.B\c\ ^ 0- Hence 

\Re{tr[paH A:Blc })\ < \tr A [pABcH A :B\c]\iM\ (10) 
= I(A:B\C)\\a\\. (11) 

It is tempting to think that the same proof holds when- 
ever all the reduced density matrices commute with each 
other. The statement itself is correct, but the preceding 
proof cannot be used. The reason why the proof worked 
for a classical systems was because there is a simulta- 
neous local product basis for all the density matrices. 
It is possible for the density operators to commute with 
each other, yet do not allow a local product basis. For in- 
stance, one can easily prove that the reduced density ma- 
trix of 2D toric code commute with each other. However, 
the simultaneous eigenstates of these operators have en- 
tanglement in general between two different subsystems. 
Proof that includes such cases is presented in Appendix 
[Bl In fact, these are the very models we are interested in. 
We denote these models as strongly commuting models. 

Definition 3. H is strongly commuting if 



[pa,Pb] = 



for all A, B. p= 

Applying the inequality from Appendix [ 
rive the following. 



(12) 



one can de- 



Lemma 2. Suppose p is a Gibbs state of strongly com- 
muting model. If nontrivial support of a is contained in 
AC or BC, 



\Re(C(H A ..B\c,v))\ < I (A: B\C)\\a\ 



(13) 
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III. MAIN RESULT 

Consider a perturbation V = J2i v * anc ^ a strongly 
commuting hamiltonian H = y ■ hj with \\hi\\, \\vi\\ < J. 
We assume UjS and hiS are sufficiently local. One way to 
state this is to use interaction radius. We say interaction 
radius of Vi to be r if its nontrivial support is contained a 
ball with radius r centered at position i. Wc shall assume 
a hnite radius of interaction r B and ry for the original 
hamiltonian and the perturbation. Main result can be 
written as the following. 



• ••/*••• 
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A 



B 



(a) 



(b) 



of trace, one can show that the term in EqfT5]is the only 
surviving one. Using Eq(51 directional derivative can be 
expressed as a connected correlation function between 
entanglement spectrum of various subsystems and local 
perturbation. 

For each of these local terms, we employ the follow- 
ing strategy. If the local term is sufficiently far away 
from the subsystems involved in computing the cntropic 
order parameter, connected correlation function can be 
bounded by C{0 A ,0 B ) < yj2I{A : B) \\0 A \\ \\ B || , where 
|| • • • || denotes an operator norm. The operator norm of 
the linear combination of entanglement entropy can be 
bounded by 0((3JVol). Operator norm of the local oper- 
ator under transform On — > Q.. ta ° x v j s nonincreasing. 
Both statements about operator norm are nontrivial. See 
Appendix [S] for detailed explanation. I(A : B) can be 
bounded by a function which depends on the distance 
between two subsystems. 

If the local term is not sufficiently far away, one can 
always deform the linear combination of entanglement 
spectrum to be at least distance l\ away from the local 
term. Connected correlation function between the de- 
formed order parameter and local term can be bounded 
similarly. An example of such deformation is depicted in 
FIGE 
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(d) 



{a)H A . D \ BC \ D (b)H A[E \ BC \( DE ) (c)-H"a:E|C\b 

FIG. 9: These three terms appear in a difference between the 
original and deformed order parameter. 

The difference between I (A : B\C) to I(A : B'\C), 
is small in a thermodynamic limit. However, we still 
need to check if the correlation between the local ob- 
servable and H a .b\c — H A -b>c is small. Terms arising 
from this difference is depicted in FIGE] Depending on 
the position of the local term, there are two possibilities 
to consider. First, the support of the local observable 
is naturally included in the support of ABC, where the 
subsystems are used for computing I (A : B\C). Since 
the distance between A and B is assumed to be at least 
I2; for large enough I2, the support of the local term ei- 
ther lies on AC or BC. Hence the correlation can be 
bounded using Lemma O It is also possible for the lo- 
cal term to have support outside of ABC. This is what 
happens in FIG El Without loss of generality, assume 
the local operator is closer to A than B. Enlarge A 
to AD such that A and D are connected. Using the 
chain rule H A:B \c = H AD , B \ C - H D:C \ A b, one can ap- 
Campbell-Hausdorff formula. Using the cyclic property ply Lemma O Such correlation function can be bounded 



FIG. 8: Effect of the perturbation can be expressed as a sum 
of connected correlation function between local terms and 
H A .B\c- h is always possible to deform the order parame- 
ter away from any local operator. 



Theorem 1. // S({Ai}) is (li^^-deformable for a 
strongly commuting hamiltonian, 

\dlS{{Ai})\ < 0(max(f(h - r ), /(| - r a )){pjVolf), 

(14) 

where r Q = r^ + ry . 

Proof. One may start with the following formula. 

d;S A = tr[-a\n PA ], (15) 
for tr[a] = 0. Simple way to see this is to use Baker- 
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Extra Assumption Improvement 
Superpolynomially decaying / Vol 2 — ► Vol 
No thermal phase transition Vol 2 — > Vol x Vol(S) 
Two assumptions combined Uo/ 2 — > Vol(S) 

TABLE I: A list of possible improvements on the bound de- 
pending on extra assumptions. The second and third means 
a combined effect of the superpolynomial decay of / and ab- 
sence of thermal phase transition. 

by 0{/3 2 J 2 f( l -f - r )Vol). There arc 0{Vol) such local 
terms. Combining all of these together, we arrive at the 
desired bound. □ 

As stated in the proof, our result applies to strongly 
commuting models. Interesting nontrivial example in 
this category is so called stabilizer models. For the proof 
of strong commuting condition, see Appendix [Cj 

There are two length scales here. First length scale 
is the size of the system, which is encoded in the vol- 
ume term of the bound. Second length scale is the size 
of the subsystems which are relevant to the computation 
of entropic order parameter. Hence one must scale the 
size of the subsystems appropriately to ensure the con- 
vergence of the bound. If / does not decay fast enough, 
such convergence may not be achieved. However, we are 
considering systems with finite correlation length. There- 
fore we expect / to decay exponentially fast, or at least 
superpolynomially. 

Under such assumption, it is actually possible to re- 
duce Vol 2 term to Vol. Note that the effect of local 
perturbations which are far away from the support of 
the order parameter was simply bounded by a product 
of the number of local terms and the maximal correlation. 
Instead of this naive bound, one can imagine taking an 
integral over this space, so that the terms that are far- 
ther away can be bounded by a small number. Additional 
factor of Vol appears from the bound \\Ha\\ < 2/3||iJ||, 
which is derived in AppendixfA] If one assume an asymp- 
totic conditional independence for all finite temperature, 
this term can be improved to 0(Vol(S)), where S is a 
union of the subsystems used for computing the entropic 
order parameter. These potential improvements are sum- 
marized in Tabic U 



IV. DISCUSSION 

One can easily see that the proof of the main result is 
based on Lemma [1] and Lemma [2] This naturally leads 
us to the question of generalizing the analogous results 
when dropping the strong commuting condition. It is 
obvious Lemma [T] will not hold for general local hamilto- 
nian models, but one can still ask if there exists a Lieb- 
Robinson type locality bound for such time evolution. 
Hastings showed such result for ID system, but a similar 
statement can be established for any hamiltonian which 



allows nontrivial Licb-Robinson bound. 

Lemma [5] is really a statement about the spectrum of 
an operator tr A\pH a-.b\c + Ha-.b\cp\ in disguise. Numer- 
ical simulation shows that this operator may have nega- 
tive eigenvalues, and existence of these negative eigenval- 
ues forbids Lemma [5] to be generalized to arbitrary den- 
sity operator. If one can bound the sum of these negative 
eigenvalues, it might be possible to establish an analogous 
result. We do not expect Lemma [2] to hold without the 
strong commuting condition, but a much looser bound 
is still applicable to the stability proof. For instance, 
poly (hid) factor, where d is the dimension of the system, 
does not ruin the proof as long as the function / decays 
sufficiently fast. 

For models which satisfy the strong commuting con- 
dition, our result resolves some of the confusing aspects 
of finite temperature topological entanglement entropy. 
At zero temperature, there are two well known ways of 
extracting topological entanglement entropy in 2D. Sur- 
prisingly at finite temperature, they show different be- 
havior. Levin- Wen configuration produces a trivial topo- 
logical entanglement entropy for a sufficiently large sys- 
tem size^S On the other hand, Kitaev-Preskill configu- 
ration produces a nonzero topological entanglement en- 
tropy that changes in temperature^ This seems puz- 
zling; two different order parameters with same predic- 
tion at zero temperature give rise to different answers at 
finite temperature. 

In the context of deformable entropic order parame- 
ter, this can be understood in a following way. First, 
Levin- Wen construction give rise to a (li, ^-deformable 
entropic order parameter for sufficiently large l\ and li. 
The same statement does not hold for Kitacv-Prcskill 
configuration. The main reason is the existence of triple 
point where subsystems A, B, and C meet. To see this, 
let us briefly review the stability proof. Finite tempera- 
ture topological entanglement entropy of Kitaev-Preskill 
configuration depends on the temperature, so let us con- 
sider a temperature perturbation (3 — > (3 + 5(3. This cor- 
responds to simply perturbing the coupling constant of 
the entire hamiltonian. 

For the local perturbation terms that are sufficiently 
far away from the triple point, it is possible to deform 
the subsystems sufficiently far away while the deforma- 
tion can be bounded by a small number. Hence the topo- 
logical entanglement entropy will be stable against local 
terms far away from the triple point. However, one can- 
not find such deformation for the terms near the triple 
point. Assuming finite correlation length, there might be 
a 0(1) term that may not be bounded unless there is an 
extra insight to get rid of it. We suspect this is the origin 
of O(l) term for 2D toric code at finite tcmpcraturci^i 
The same logic can be applied generally. 

This is an evidence that the configurations that can 
extract topological entanglement entropy at zero temper- 
ature may not be a suitable choice at finite temperature. 
Assuming such cancelation occurs for a set of subsys- 
tems {^4i}, these phenomena occur whenever there is a 
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triple point that at least three of these subsystems meet 
together. For this reason, we expect stacking Kitaev- 
Prcskill type construction will not be a meaningful or- 
der parameter that can detect topological phase. Higher 
dimensional generalization of Levin- Wen construction 
seems to avoid this problem, as evidenced by the exam- 
ples presented in this paper. 



V. CONCLUSION 

We presented a general formalism to prove a stability of 
finite temperature topological entanglement entropy un- 
der hamiltonian perturbation in a general context. When 
applied to most of the exactly solvable models exhibit- 
ing topological order at finite temperature, it shows that 
the topological entanglement entropy is invariant under a 
sum of geometrically local bounded- norm terms in a ther- 
modynamic limit. We also presented a general reason to 
believe why certain topological entanglement entropy at 
zero temperature may not be stable against hamiltonian 
perturbation at finite temperature. 

The technical foundation of this work is based on two 
statements. First statement concerns a certain variant of 
Licb-Robinson type bound in quantum many body sys- 
tems. Second statement concerns the spectrum of some 
operator whose trace is conditional mutual information. 
The first statement in fact holds generally for any quan- 
tum many-body system which satisfies Licb-Robinson 
bound. Generalization of second statement seems elusive 
at this point. This is the only missing link to proving the 
stability of topologically invariant order parameters at fi- 
nite temperature. Whether such statement is true or not 
remains as an open problem. 

Another interesting question concerns the condition we 
imposed on the state, which is the asymptotic conditional 
independence. Circuit definition of topological order in 
rough terms, assumes a local indistinguishability of topo- 
logically different sectors. It is clear that the asymptotic 
conditional independence does not imply the circuit def- 
inition of topological order, since the 2D finite tempera- 
ture examples presented in this paper are not topologi- 
cally ordered. It would be interesting to see if the circuit 
definition of topological order implies asymptotic condi- 
tional independence. 
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Appendix A: Useful facts from Matrix Algebra 

Let H n be a n-dimensional Hilbcrt space. A set of 
operators acting on this Hilbert space is denoted as 
B(H n ) : "Hn — > Hn- We are interested in a specific set 
of superoperators S(H n ) : B(H n ) — > B(H n ) that share a 
common property. 

A family of superoperators we study are the ones char- 
acterized by a function / : R d — > C, /(0) = 1 that maps 
an operator O in a following way. 

(),, -O.Jl.r, .r,i (Al) 

where i,j corresponds to the row and column of the ma- 
trix and Ei G (— oo, oo). The main question is to find the 
necessary and sufficient condition for such map to be a 
quantum operation. One can easily see that the trace is 
preserved by the the property /(0) = 1. It is tempting to 
think that such maps are quantum operation, since the 
entries of the matrices contract under the map. However, 
this is not necessarily true. For instance, consider a map 
with f(xi — Xj) = f{xj — Xi) = and f(x) — 1 every- 
where else. One can easily check that it maps a matrix 
Oij = 1 Vi,j = 1,2,3, which has eigenvalue 3,0,0 to 
another matrix that has negative eigenvalue. 

A systematic way of deciding if the map is a quantum 
operation or not, is to use Choi's theorem on completely 
positive maps. Applied to these superoperators, Choi's 
theorem tells us that our map is completely positive if 
and only if the following matrix is positive semi-definite. 

F ij =f(x i -x j ) (A2) 

It was Bochner who discovered a necessary and and 
sufficient condition for the positivity of such matrices. 
In our setting, Bochner's theorem implies the following 
statement^ 

Theorem 2. (Bochner)Let f : M. d —> C be a continuous 
function. F^j is positive semi-definite if and only if there 
exists a positive function g such that 

f(x) = J e 2 ^ x g{p)d d p. (A3) 

This is the standard formulation of Bochner's theorem, 
but the same line of standard proof holds even when g(p) 
is a distribution as long as it is positive^ Easiest exam- 
ple to see this is the unitary operation. The correspond- 
ing function is f{x) = e lx . The inverse Fourier transform 
of this function is a delta function. 

Let us carry on with another example. Consider 
f{x) = sin i^ x j ■ The inverse fourier transform is g(p) = 

^3/2 sech 2,7rps 

2 ^ 2 . This may seem like a contrived example, 

but a quantum operation associated with this function 
naturally arises when one studies the entanglement spec- 
trum of a finite temperature equilibrium statc(Gibbs 
state). Suppose we have a hamiltonian H. Gibbs state at 
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finite temperature is p = ^-g-, where Z = tr(e~ H ). Typ- 
ically one writes (3H in the place of H, where (3 being 
the inverse temperature. However, we shall just absorb 
it into a definition of H for notational convenience. Wc 
would like to study the effect of perturbation H — > H+eV 
on the entanglement spectrum. One can easily check di- 
rectional derivative of p A can be computed in a following 
way. 

-B%tr A e [e- H ] = tr A , [ f e^Ve^ 1 "^] (A4) 
Jo 



tr A ,[e^ H / e- sH Ve sH e-i H ds] 



tr A a [e-5 H $(F)e _ 4 



(A5) 
(A6) 



where we used Duhamel's formula on the first line. From 
the first line to the second line we simply rearranged the 

terms. On the third line, $(V) = J\ e~ sH Ve sH ds. If 

2 

we expand V in terms of the eigestates of H, $ maps 
Vij —> Vij smi ^- x ) ; where x = ^'"^ and Ei is the zth 
energy eigenstate. In other words, $ is an inverse of a 
quantum operation <&(V)ij = Vij sin ^ x ^ ■ As one can sec 
from this example, one can define a set of superoperators 
defined by a hamiltonian H and a function f(x). Wc 
shall denote such operation and its inverse as following. 



Definition 4. 



•if , : On o,, 



(A7) 
(A8) 



Oij = (i\0\j), where \i) is the ith eigenstate of H and 
Ei is the corresponding eigenvalue. 

From these results, we can show the following lemma. 

Lemma 3. $f( x \ is a quantum operation if and only if 
f(x) has a positive inverse Fourier transform. 

Corollary 1. 

\\<n^/*(V)\\ < \\V\\, (A9) 
Proof. It follows from following two identities. 



tanha; 



x 



E 

fc=0 



X 2 + (k+ ±) 2 7T 2 



(A10) 



Lemma 4. For \\H\\, \\V\\ < oo, 

d p H A = (H) - $f A (H) A a {A , (A12) 

where f(x) = sin ^/ 2 ) ■ (°)a°\a = P A *tr A c \p^Op^\p A ~- 
is a quantum conditional expectation. 



Corollary 2. 



||tfA||<2/?pr| 



(A13) 



Proof. $f A is a quantum operation. Conditional expec- 
tation is also a quantum operation^ Using the norm 
contractivity of quantum operation, one can easily ob- 
tain \\H A \\ < 2/?||F||. □ 

In the absence of finite temperature phase transition, 
one can expect an existence of finite length scale at all 
finite temperature. In this case, it is possible to cancel 
out most of the correlations residing far away from A. 
We can get a better bound this way, but for the purpose 
of the paper this is not so crucial. 



Appendix B: Modification of Klein's inequality 

Klein's inequality states the following statement. 

Theorem 3. For A, B > , 

tr[AQnA-]nB)]>tr[A-B]. (Bl) 

What wc would like to prove is an analogous statement 
for partial trace. 

Lemma 5. For A, B > 0, [A, B] = 0, 

trc [A(\n A -In B)\ > tr c [A - B] (B2) 

Proof. Let \i) denote a common eigenstate of A and B. 
a% and bi denote the corresponding cigenstates. pi = 
tr c [\i) 

tr c [A{\nA-\nB)\ =^a i (ln^)p i (B3) 



>]>>(1-^V 4 (B4) 

* * ft. A 



tr c [A-B}. 



(B5) 



From the first two the second line, we used In - > 1 — 
x. □ 



1 + x 2 V 2 



(All) 
□ 



Corollary 3. For the Gibbs state of strongly commuting 
hamiltonian, 

tr A [p AB c{H A BC + H B - H BC - H AB )] > 0. (B6) 
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Appendix C: Strongly Commuting Models 

Strongly commuting condition seems to be a very 
strong constraint, but fortunately there are number of 
interesting models that satisfy this condition. First class 
of models we would like to consider are models based on 
quantum error correcting code. Stabilizer code S is an 
abelian subgroup of a generalized pauli group V n that 
does not contain —I as a group element. If there exist a 
spatially local group generators for S, one can associate it 
to a quantum many-body system with hamiltonian corre- 
sponding to those generators. Since the group is abelian, 
the local generators commute with each other. Further- 
more, all the elements of the group is hermitian, so the 
corresponding hamiltonian is simply a sum of the gener- 
ators. Furthermore, since the stabilizer group elements 
are vanishes under any partial trace over its nontrivial 
support. Using this fact, one can prove the following. 



Lemma 6. pa = Es es(yi) c i^' for some coefficients 

{a}- 

Proof, p can be expanded as a sum of stabilizer group 
elements. When taking the partial trace, any operator 
that has nontrivial support on A c vanishes. The ones 
that have nontrivial support only on A survives, and they 
are again a stabilizer group element. Hence pa can be 
written as a sum of elements in S(A). □ 

Corollary 4. For Gibbs state of stabilizer hamiltonian, 

[pa,Pb] = 0. (CI) 

Stabilizer models cover interesting models in many dif- 
ferent dimensions, including Toric code in 2,3,4 spatial 
dimension, topological color code, quantum glass code, 
and other variants of 3D models. 
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